Charged-Particle Pseudorapidity Distributions in Au+Au Collisions at y/s^N = 62.4 
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The charged-particle pseudorapidity density for Au+Au collisions at y / 5jvjv=62.4 GeV has been 
measured over a wide range of impact parameters and compared to results obtained at other energies. 
As a function of collision energy, the pseudorapidity distribution grows systematically both in height 
and width. The mid-rapidity density is found to grow approximately logarithmically between AGS 
energies and the top RHIC energy. As a function of centrality, there is an approximate factorization 
of the centrality dependence of the mid-rapidity yields and the overall multiplicity scale. The new 
results at y / ijvjv=62.4 GeV confirm the previously observed phenomenon of "extended longitudinal 
scaling" in the pseudorapidity distributions when viewed in the rest frame of one of the colliding 
nuclei. It is also found that the evolution of the shape of the distribution with centrality is energy 
independent, when viewed in this reference frame. As a function of centrality, the total charged 
particle multiplicity scales linearly with the number of participant pairs as it was observed at other 
energies. 
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In previous publications the PHOBOS Collaboration 
has presented the full systematic behavior of inclusive 
charged particle production in heavy ion collisions over a 
large range of collision parameters: 1) ^snn from 19.6 
GeV to 200 GeV, 2) Pseudorapidities from r] = -5.4 to 
5.4, nearly the full solid angle, and 3) average impact pa- 
rameter from (b) = 3 — 10.5 fm, corresponding to 50-360 
participating nucleons {N part ) QEIJSEIiH- While 
these data are already useful as a broad systematic study, 
several non-trivial features have been made manifest by 
direct comparisons between the data at different energies 
and centralities: 1) a logarithmic increase with ^/snn in 
the mid-rapidity particle density 0, 2) an approximate 
factorization of the centrality and energy dependence of 
the mid-rapidity yields Q 3) the phenomenon of "limit- 
ing fragmentation" in the forward direction [fjj , and 4) a 
linear " V part -scaling" of the total particle yield || . This 
paper presents for the first time the multiplicity data 
for Au+Au collisions at the most recent RHIC energy of 
y/SNN = 62.4 GeV, corresponding to the top energy re- 
ported by several earlier p + p experiments at the CERN 
ISR. With our new data, we can test the previously-found 



scaling relationships at an intermediate energy. 

The data were obtained with the PHOBOS detector0 
during the RHIC 2004 run. In this analysis only the 
data taken with the magnetic field switched off are an- 
alyzed. To select events with various ranges of impact 
parameter, characterized typically by N part , we use the 
particle multiplicity measured in two sets of "paddle" 
counters, situated at z = ±3.21 meters from the nominal 
interaction point. These cover a pseudorapidity range of 
3.2 < |»7| < 4.5 with 95% azimuthal acceptance. The 
overall triggering and event selection efficiency at 62.4 
GeV corresponds to 81 + 2% of the total Au+Au inelastic 
cross section, estimated using HIJING simulations 10]. 
We use the Glauber model calculation implemented in 
HIJING to estimate (N part ) for each centrality bin by 
assuming a monotonic relationship between N part and 
the relevant experimental observable. This procedure has 
been described in Refs. 0,0,0- It was found that trigger 
efficiencies are typically around 100% for the top 50% of 
the total cross section at energies of ^snn = 62.4, 130 
and 200 GeV. At the lowest energy of ^snn — 19.6 
GeV, an alternative method of centrality determination 
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FIG. 1: Pseudorapidity distributions dN c h/dr) from Au+Au 
collisions at ^/snn = 62.4 GeV. Open symbols show the re- 
sults obtained using the paddle-based centrality method with 
90% C.L. systematic errors indicated by grey bands for se- 
lected bins. The dotted lines show the results obtained using 
the Octagon-based centrality method. The filled circles show 
the average value of dN c h/ 'd7/j|, ; | <:L using the tracklet tech- 
nique, with the horizontal bars indicating the interval in r\ 
over which the averaging is done. The centrality is denoted 
by the fraction of the total inelastic cross section, with smaller 
numbers corresponding to more central events. 

was developed, that uses all particles detected in the "Oc- 
tagon" silicon detector, covering |7y| < 3, as a measure of 
the particle multiplicity Jfj. These two methods agree 
within 3% at the higher energies and have been used as a 
cross-check in this analysis. However it should be noted 
that the two methods give values of (N par t) that differ 
by 2% in the most central bin, since this bin is more sen- 
sitive to the fluctuations of the variable used to estimate 
the centrality. Thus, we separate the tabulated results 
(in Table HJ) for that bin, but do not in the other bins 
when they agree to better than 1%. 

Several methods were used to estimate the charged 
particle density in each centrality class. In the full phase 
space, a combination of data from the single-layer "Oc- 
tagon" (|ry| < 3) and "Ring" (3 < \r)\ < 5.4) detectors 
were analyzed using two different techniques. In the 
"analog" method, the deposited energy in a detector pad 
is used to estimate the number of particles traversing the 
pad after accounting for orientation of the Si-wafer rel- 
ative to the interaction point. The "digital" approach 
treats each pad as a binary counter and assumes Pois- 
son statistics to estimate the total occupancy in vari- 
ous regions of pseudorapidity. These methods have been 
discussed in more detail in Refs. Q and Q. At mid- 
rapidity, the PHOBOS vertex detector, consisting of two 
planes covering \rj\ < 0.92 over a limited azimuthal range, 
A(f) ~ 90°, is used to count "tracklets". These are two- 
hit tracks which point back to the event vertex, providing 
redundancy not present in the single-layer analyses, and 
thereby reducing systematic effects at mid-rapidity. This 



FIG. 2: The pseudorapidity distributions measured in the 6% 
most-central Au+Au collisions at four RHIC energies. 90% 
C.L. systematic errors are shown as grey bands. 

method has been described in detail in Refs. 0,0,0- 

The y/sNN = 62.4 GeV data are shown in Fig. ^ as a 
function of collision centrality, determined by the paddle- 
counter method. Results using the Octagon-based cen- 
trality method, shown for each bin by dotted lines, agree 
very well with those from the paddle-based method. Mid- 
rapidity data from the tracklet method also agree well 
with the single-Si-layer analysis over the full centrality 
range studied. 

To place these data in context, Fig.[2]shows data from 
the 6% most central collisions in comparison with sim- 
ilar data at 19.6, 130, and 200 GeV from Ref. @. In- 
creases both in the height and width of the distribution 
are observed as a function of increasing energy. Already 
at ^/sjvat = 62.4 GeV, the central "plateau" nascent 
at ^/snn — 19.6 GeV is fully developed and grows in 
width slowly with increasing energy. Of course, the ex- 
istence of a plateau in dN c h/ drj does not necessarily im- 
ply the existence of a plateau in dN/dy because of the 
non-trivial transformation between rapidity and pseudo- 
rapidity (dy = fidrf) [6j. 

The scaled particle densities near midrapidity 
{d^Jd^^lN^riM) for = 19-6, 130 and 200 

GeV [J H H H HQ, El are shown as a function of col- 
lision energy in Fig. |3| for the 6% most central events. 
Where possible, PHOBOS results from the various mea- 
surement techniques have been averaged at each energy, 
weighted by the inverse square of the relative error. At 
^snn = 62.4 GeV, this gives dN ch /di]\ M<1 / (N part /2) = 
2.64+ 0.18. Data from comparable centralities at lower 
energies are also shown, as compiled in Ref. H3. While 
high-statistics data from the three other RHIC energies 
suggested an approximately-logarithmic rise of the par- 
ticle density, the low-statistics data point measured at 
t/snn — 56 GeV was found to be only barely consis- 
tent with a logarithmic fit based only on the data at 
•J snn = 19.6, 130 and 200 GeV (shown as a dashed line). 
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FIG. 3: dN c h/dr]\\ rt \ < -i/{N P art/'2) shown for Au+Au collisions 
as a function of energy. The PHOBOS data, averaged over 
all available measurement techniques, is compared with lower- 
energy A+A data as well as a variety of p + p and p + p data. 
The thick dashed line is a fit (a + 61n(s), with a = —0.40 and 
b = 0.39) to the ^^sWn = 19.6, 130 and 200 GeV data points. 
The inelastic p + p and p + p data have been fit by a function 
a + bs c (with a = 0.35, b = 0.52 and c = 0.12), shown by a 
thin dotted line, used only for interpolation. 
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The data point at ^/snn — 62.4 GeV (indicated by an 
arrow in Fig.[3J) falls closer to the fit, and is, within errors, 
consistent with the logarithmic rise. The A + A data are 
compared to a wide range of p + p and p+p data, sepa- 
rately shown for inelastic as well as non-single diffractive 
events ^2 ^| . While these data also appear to rise 
logarithmically at higher energies, the inclusion of data 
points below y/s ~ 30 GeV appears to indicate a curva- 
ture in the energy dependence. To interpolate between 
the measured points, they have been fit by a function 
a + bs c , shown by the thin dashed line. 

The centrality dependence oidNch/drjl^i^/ (N part /2), 
measured with the tracklet technique, is shown in Fig.^i, 
and tabulated in Table [i] Data at ^/snn = 19.6,62.4, 
and 200 GeV have been analyzed using only the track- 
let technique in the PHOBOS Vertex detector (hereafter 
called the "vertex-tracklet" method) and the Octagon- 
based centrality method. At 130 and 200 GeV, results are 
also available using a method (the "combined" method) 
which averages the vertex tracklet results with a similar 
method using the PHOBOS spectrometer and us- 
ing the paddle-based centrality method. The results are 
typically compatible within 2% over the full centrality 
range, as can be seen by direct comparison in Fig. 0^ of 
the vertex tracklet and the combined result for 200 GeV. 
The vertex-tracklet method and Octagon-based central- 
ity method is used at 62.4 GeV for overall consistency 
and partial cancellation of certain systematics in the ra- 
tios relative to 19.6 and 200 GeV data. 

The centrality dependence of the mid-rapidity yields 
has often been interpreted in a two-component picture 



FIG. 4: a.) dN c h/dT]\\ ri i < i/(Npart/2) measured using the 
tracklet technique as a function of N par t for four RHIC ener- 
gies (y/s^ = 19.6, 62.4, 130 and 200 GeV). Also shown are 
inelastic p + p data for y/s = 62.4 and i/s — 200 GeV, and in- 
terpolated values for = 19.6 GeV (Ref. 0) and = 130 
GeV (using the fit shown in Fig. |3J. The systematic errors are 
shown as 90% C.L. ellipses (reflecting the trivial correlation 
in dN ch /dri\^i <1 /(Npart/2) and N part ). b.) The ratio of the 
scaled pseudorapidity densities measured at 200 GeV to those 
measured at ^/sjvjv = 19.6, 62.4 and 130 GeV. 



of particle production, with "soft" processes scaling with 
Np ar t and "hard" processes scaling with the number of 
binary collisions, N co ll- As the beam energy increases, 
particle production from the hard processes, which ex- 
ceed the number of participant pairs by a factor of ~ 5 — 6 
in central events for ^snn ranging from 19.6 to 200 GeV, 
are expected to dominate over that from soft processes 
as the minijet cross sections increase [lOj. This expec- 
tation may be examined by studying the ratio of the 
yields at different energies for the same fraction of the 
total cross section, shown in Fig. Despite the ex- 
pected increase in hard processes with increasing energy, 
these ratios are observed to be constant over the mea- 
sured centrality range, showing a "factorization" of beam 
energy and collision geometry at midrapidity. This re- 
sult extends the analysis presented in Ref. and is fully 
compatible with the constant ratio found in comparisons 
between data at ^/snn — 200 and 130 GeV obtained 
with the "combined" method. All of these data suggest 
that while two-component models can fit the midrapidity 
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FIG. 5: The scaled pseudorapidity density 

(dN c h/dri/{Npart/2)) for two centrality bins (0-6% and 
35-40%) and four RHIC energies shown in the rest frame of 
one of the projectiles by using the variable rj' = rj — ybeam- 
For clarity, systematic errors are not shown in this figure. 



data at each energy, they do not have an energy depen- 
dence characteristic of a growing contribution of hard- 
processes. 

In a previous PHOBOS publication, the phenomenon 
of "limiting fragmentation" was observed by comparing 
Au+Au collisions at the three RHIC energies for which 
An data were available This phenomenon, which we 
also refer to as "extended longitudinal scaling" is sim- 
ply the invariance of the scaled yields dN c h/drj/ (N part /2) 
with beam energy in the reference frame of one of the pro- 
jectiles, i.e. by plotting the scaled yields with respect to 
the variable rf = r\ — yteam, where ybeam is the beam 
rapidity. The concept of extended longitudinal scaling 
is expected to apply to rapidity distributions, but since 
q r~i y — ln(pT /tut) for particles emitted far away from 
mid-rapidity, this scaling is also expected to apply to 
pseudorapidity distributions. This scaling phenomenon 
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0-3% 


356±11 


492±36 


2.76±0.23 








349±11 






2988±149 


17.10±1.02 


3-6% 


325±10 


433±32 


2.67±0.22 


2775±138 


17.17±1.01 


6-10% 


288± 9 


377±28 


2.62±0.21 


2489±124 


17.29±1.01 


10-15% 


248± 8 


316±23 


2.55±0.21 


2120±106 


17.13±1.00 


15-20% 


209± 7 


260±19 


2.50±0.20 


1777±88 


17.03±1.02 


20-25% 


174± 7 


212±15 


2.44±0.21 


1485±74 


17.07±1.08 


25-30% 


145± 7 


174±13 


2.41±0.21 


1236±61 


17.03±1.17 


30-35% 


119± 7 


140±10 


2.35±0.22 


1027±51 


17.15±1.30 


35-40% 


98± 7 


111± 8 


2.28±0.23 


840±42 


17.17±1.45 


40-45% 


78± 6 


87± 6 


2.24±0.25 


679±33 


17.30±1.63 


45-50% 


62± 6 


67± 5 


2.16±0.26 


532±26 


17.16±1.80 



TABLE I: Data for Au+Au collisions at y?^v = 62.4 GeV 
including dN ch /dr]\ M<1 , dN ch /dr]\ M<1 /{N part /2) as shown 
in Fig. 4 and (N ch ) and (N ch }/(N par t/2} as shown in Fig. 7. 
The difference in N par t in the — 3% bin is explained in the 
text. 



was also found in d+Au collisions at RHIC 0] and, sur- 
prisingly, also for elliptic flow at all of the RHIC energies, 
including ^/snn = 62.4 GeV The new multiplicity 
data at ^snn = 62.4 GeV, shown in Figs. andEt, 
clearly fit well into the existing pattern. 

The data in Fig.[5]suggest that the change in the shape 
of dN c h/d.T] away from mid-rapidity is a stronger function 
of the collision geometry than of the beam energy, when 
observed in the rest frame of one of the projectiles. This 
is shown more clearly in Fig. El where the ratio of pe- 
ripheral to central data scaled by N part 



i?p?T(7/,35 - 40%) 



FIG. 6: Rp p c art as a function of rf comparing the 35-40% bin 
to the 0-6% most central bin for four different energies. 



is plotted as a function of rf . The error bars in this fig- 
ure indicate 90% C.L. systematic errors. A similar figure 
was shown in Ref. 0, showing that the change in shape 
as a function of centrality is independent of beam energy 
when observed in the rest frame of one of the projectiles. 
This ratio was also used to compare the shapes of trans- 
verse momentum distributions, dN/dpx, measured near 
mid-rapidity for different energies and centralities, and 
a similar invariance with energy was found for each cen- 
trality bin [rj|. i.e. in all of these cases, the centrality 
and energy dependences factorize to a large extent. Such 
behavior has also been seen in proton-nucleus collisions 
at lower energies |T^ . 

As observed previously, the centrality dependence 
of the limiting curve has the interesting property 
that the decrease in the scaled particle density, 
dN c f l /drj/(N par t/2} 1 at mid-rapidity when moving from 
central to peripheral events is correlated with the in- 



■5 



CM 



£30 



20 



10 




200 GeV 



130 GeV 



4 GeV 



= 4 $ £ £ $£19.6 GeV 





P(P)+P 



200 400 
N 



part 



FIG. 7: (N c h) I (Npart /2) , obtained by extrapolating the data 
at each energy into the unmeasured region, as a function of 
centrality. The 90% C.L. uncertainty on N c h and N par t have 
been combined into the error bars, while the 90% C.L. uncer- 
tainty on the extrapolation procedure is indicated by a grey 
band. Inelastic p + p and p + p data, interpolated using a 
power-law fit, are shown at N par t = 2. 



crease at forward rapidities. Although some of the par- 
ticles with rf > may be attributed to emission from 
spectators, the systematic change of the slope leads to 
an approximately constant total multiplicity. Using the 
method outlined in Ref. 8] , which combines analytic fits 
of the measured region with estimates of the unmeasured 
yields using the lower-energy data, the total charged- 
particle multiplicity extrapolated to 47r has been calcu- 
lated as a function of centrality, as shown in Fig. {7\ and 
tabulated in Tabled The uncertainty from the extrapo- 
lation procedure itself is indicated by the grey bands. As 
at the other RHIC energies 0, the data at ^/snn — 62.4 
GeV shows an approximately linear relationship between 
N c h and Npart. ■ This persistence of "wounded nucleon 
scaling" |l9|, |20j has not been fully explained for heavy 
ion collisions, especially since the multiplicity is evidently 
not a simple multiplication of nucleon-nucleon multiplic- 
ity by N part /2. 

It is an interesting question whether or not the vari- 
ous scaling behaviors discussed, such as the factorization 
of energy and geometry at midrapidity (Fig. and the 
similar factorization of the distributions in rj' (Fig. EJ, 
should be considered independent phenomena. Already, 
the N part scaling shown in Fig. suggests that modifi- 



cations to particle production at forward rapidities are 
strongly correlated with compensating changes at mid- 
rapidity. If, in fact, the pseudorapidity distribution at 
each energy deviates from the limiting curve at around 
rj r*> 1.5 — 2 by flattening out, and this deviation is cen- 
trality independent (a statement which is broadly consis- 
tent with the available data, except perhaps at the lowest 
RHIC energy) , then the factorization of energy and geom- 
etry at mid-rapidity follows naturally as a consequence of 
the centrality dependence of the energy-independent lim- 
iting curve. The same condition also relates the approx- 
imately logarithmic energy dependence (Fig. |3J to the 
shape of the limiting curve in central events. Of course, 
empirical observations like this do not explain why these 
relationships between the various regions of phase space 
hold, but rather point to issues that need to be addressed 
in the global understanding of heavy ion collisions. 

In summary, the charged-particle pseudorapidity den- 
sity has been measured by PHOBOS for Au+Au colli- 
sions at y/sNN~ 62.4 GeV, matching the top ISR en- 
ergy. As a function of collision energy, the pseudorapid- 
ity distribution grows systematically both in height and 
width. The mid-rapidity density is found to grow approx- 
imately logarithmically between AGS energies and the 
top RHIC energy. As a function of centrality, there is an 
approximate factorization of the centrality dependence of 
the mid-rapidity yields and the energy-dependent overall 
multiplicity scale. The phenomenon of "extended lon- 
gitudinal scaling" (also known as "limiting fragmenta- 
tion") is clearly present in the 62.4 GeV data. Finally, 
a relatively-small extrapolation of the measured yields 
to 47T allows the extraction of the total charged-particle 
multiplicity. As at the other RHIC energies, N c h is found 
to scale approximately linearly with the number of par- 
ticipants over the range of collision centralities studied. 
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